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In a comment on arXiv:1006.5070vl, Drechsler et al. present new band-structure calculations 
suggesting that the frustrated ferromagnetic spin- 1/2 chain LiCuVCU should be described by a 
strong rather than weak ferromagnetic nearest-neighbor interaction, in contradiction with their 
previous calculations. In our reply, we show that their new results are at odds with the observed 
magnetic structure, that their analysis of the static susceptibility neglects important contributions, 
and that their criticism of the spin-wave analysis of the bound-state dispersion is unfounded. We 
further show that their new exact diagonalization results reinforce our conclusion on the existence 
of a four-spinon continuum in LiCuVCU, see Enderle et al, Phys. Rev. Lett. 104 (2010) 237207. 



In their Comment, Drechsler et al. [lj argue that the 
frustrated ferromagnetic spin-1/2 chain LiCuV04 should 
be described by a strong ferromagnetic (FM) nearest 
neighbor (NN) interaction J\ , rather than a weak FM J\ 
as found by an analysis of the static susceptibility and 
the dispersion of the bound state @ and also from an 
analysis of the intensity of the continuum scattering JjJ. 
Drechsler et al. have in [l[ revised their band-structure 
calculations of Q and now propose two alternative sets 
of exchange integrals with strongly FM J\, which are 
claimed to be in [better] agreement with experiment. The 
set A has J\ = —6.3 and J 2 = 5.05 meV while set B has 
J\ = —8.8, J 2 = 6.5, and an antiferromagnetic (AFM) 
interchain coupling (cf. Q) J 4 = 0.5 meV. 

We first note that the large AFM J4 of set B leads to a 
propagation vector with a component along the a* axis, 
in contradiction with the observed magnetic structure. 

Secondly, the inverse static susceptibility of 0] was re- 
analyzed in [l[ by fitting a linear relation to the high- 
temperature part above 500 K, with the (/-factor and 
Ocw as free parameters. Figure la shows the Curie- Weiss 
temperature Ocw from such linear fits of the same data 
Q as a function of the lowest temperature T m - m used 
(triangles). Clearly, the cut-off at T m i n =500 K used by 
Drechsler et al. [l| (dashed line) is fully arbitrary and 
their Ocw is not unique. If the susceptibility is corrected 
for the known constant diamagnetic and Van-Vleck sus- 
ceptibilities and the g-factor fixed to the value precisely 
known from ESR [3] , dew is rather independent of the fit 
interval (circles) and always negative, in agreement with 
the exchange integrals of 0, Q , but not with set A or B 
of [if. Additionally, our exchange integrals, where J\ is 
weakly FM, agree with the more sophisticated analysis of 
the static susceptibility performed in [i[ , which included 



high-temperature series expansion and 7V=16-ring calcu- 
lations. 

Thirdly, Drechsler et al. further argue that it is not pos- 
sible to find a unique set of J\ and J 2 from the dispersion 
of the main peaks in inelastic neutron scattering (INS) 
data. Within spin-wave (SW) theory, this statement is 
incorrect. Figure lb shows that the SW dispersion along 
the chain direction depends strongly on the choice of J\ 
and J 2 . We note that the red curve shown in Fig. lb,c 
of [l| does not correspond to the SW description ]H . 

The exact diagonalization [l[ reveals that the intensity 
maximum at energy E m / J2 and k w 37r/4 varies little 
between Ji/J 2 = -0.45 and -1.25 (E m /J 2 « 1.67 and 
1.41). If we use this exact diagonalization result to esti- 
mate J 2 directly from the experimental value E m =6.2 
meV we find J2=3.8-4.4 meV in agreement with 0, 0]. 
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FIG. 1: (a) Curie- Weiss temperature Ocw as a function of 
the lowest temperature T m i n used for fits of the inverse sus- 
ceptibility. (The highest temperature is always 640 K). (Blue) 
triangles: Procedure of [l| . (Red) circles: Fits after correction 
for van-Vleck and diamagnetic contributions with g a fixed to 
ESR value [j. (b) Chain dispersion compared to spin- wave 
calculations with different exchange parameters. 
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J 2 as large as 5 meV or even 6.5 meV can clearly be ex- 
cluded. The INS data reveal considerable intensity above 
an upper two-spinon boundary ttJ 2 sin k calculated with 
even the largest J 2 = 4.4 meV. Interestingly, such in- 
tensity above 7rJ2sinfc is also visible in the exact diag- 
onalization results for both weak (J1/J2 = —0.45) and 
strong (J1/J2 = -1.25) FM coupling Q. This inten- 
sity is in both cases much stronger than in the AFM 
Heisenberg chain For weak FM coupling, this con- 
firms four-spinon excitations above a two-spinon contin- 
uum, and reinforces our conclusion Q on the existence 
of a four-spinon continuum in LiCuVO,}. 

In our RPA model of two coupled Heisenberg chains 
Q a significant modification of the isolated Heisenberg 
chain two-spinon continuum is found for FM coupling 
between the chains, in good agreement with the main 
experimental features. The fit of both J\ and J 2 leads 
to a quantitatively satisfactory overall description of the 
observed intensity with \J\\ < \J 2 \ and values close to 



those found in [2| . The RPA-approach appears therefore 
justified, at least a posteriori. 

In conclusion, the parameter sets suggested in the 
Comment [l[ are at odds with several experimental prop- 
erties, such as the susceptibility, the magnetic ordering 
vector for one of the parameter sets, and the observed 
dispersion of the intensity maxima 0, while the param- 
eters proposed in 0, Q lead to a consistent description of 
these properties. 
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